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Abstract 

We investigate the dynamical behavior in the large scale region of non-equilibrium sys- 
tems, by employing data on the assessed value of land in 1983 - 2006 Japan. In the system 
we find the detailed quasi-balance, which has the symmetry: x\ h a j 2 9 {xi and X2 are 
two successive land prices). By using the detailed quasi-balance and Gibrat's law, we derive 
Pareto's law with varying Pareto index annually. The parameter corresponds with the 
ratio of Pareto indices (/xi + l)/(/i2 + 1), and the relation is confirmed in the empirical data 
nicely. 
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1 Introduction 



In the large scale region of wealth, income, profits, assets, sales, the number of employees 
and etc (x), the cumulative probability distribution P(> x) obeys a power-law: 

P(> a?) oc aT** . (1) 

This power-law and the exponent \x are called Pareto's law and Pareto index, respectively [1]. 

Recently, Fujiwara et al. [2] have explained Pareto's law (and the reflection law) by using 
the law of detailed balance and Gibrat's law [3], which are also observed in empirical data. The 
detailed balance is time-reversal symmetry observed in a relatively stable economy: 

Pl2(.Xl,X 2 ) = Pu(x2,Xl) . (2) 

Here x\ and x 2 are two successive incomes, profits, assets, sales, etc. and Pi2(xi, x 2 ) is a joint 
probability distribution function (pdf). On the other hand, Gibrat's law is valid in the large 
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scale region where the conditional probability distribution of growth rate Q{R\x\) is independent 
of the initial value x\\ 

Q(R\ Xl ) = Q(R) . (3) 

Here growth rate R is defined as the ratio R = xijx\ and Q(R\x\) is defined by using the pdf 
P\(x\) and the joint pdf Pm(xi,R) as Q(R\x\) = P\r{x\, R)/P\(x\) . In the proof, Fujiwara 
et al. assume no model and only use these two underlying laws in empirical data. In Ref. [4], it 
is reported that the Pareto index is also induced from the reflection law. 

These findings are important for the progress of econophysics. Above derivations are, how- 
ever, valid only in the economic equilibrium where the detailed balance (2) holds. In order to 
discuss the transition, the dynamics should be established by investigating long-term economic 
data in which dynamical transitions are observed. Unfortunately, it is difficult to obtain personal 
income or company size data for a long period. 

In this study, we investigate the dynamical behavior in the large scale region of non- 
equilibrium systems, by employing data on the assessed value of land in 1983 - 2006 Japan. 
Because the distribution of Japanese land prices has similar features with one of personal in- 
come and company size [5], and the long-term database is readily available [6]. 

In the non-equilibrium system we find the detailed quasi-balance, which has the symmetry: 
xi «-> a x-i ■ By using the detailed quasi-balance and Gibrat's law, we derive Pareto's law with 
varying Pareto index annually. The parameter 6 corresponds with the ratio of Pareto indices 
(/Lti + l)/(/U2 + 1), and the relation is confirmed in the empirical data nicely [7]. 

2 Detailed quasi-balance 

In Japan, land is a very important asset and land prices change annually in a 24-period 
(1983 - 2006). This period contains bubble term (1986 - 1991) caused by the abnormal rise of 
land prices. The economy correlates with land prices. We employ the database of the assessed 
value of land, which indicates the standard land prices, covering the 24-year period from 1983 
to 2006. 1 

The cumulative probability distributions of land prices are shown in Fig. 1-5. From Fig. 1 - 

5, the power-law is confirmed in the large scale region. For each year, we estimate Pareto index 

[i in the range of land prices from 2 x 10 5 to 10 7 yen/m 2 where Pareto's law holds approximately. 

Annual change of Pareto index fi from 1983 to 2006 is represented in Fig. 6. In this period, 

Pareto index has changed annually. This means that the system is not in equilibrium and the 

detailed balance (2) does not hold. Actually, the breakdown in the large scale region is observed 

in the scatter plot of all pieces of land assessed in the database (Fig. 7 for instance). There is 

no x\ <-> X2 symmetry in Fig. 7 obviously. On the other hand, the detailed balance (x\ <-> X2 

symmetry) in the large scale region is observed approximately in Fig. 8 for instance. 

1 In Ref. [7], the number of data points of land prices increased gradually, because the database only contained 
data points which existed in the 2005 evaluation. In this study, the database contains all data points which existed 
in every year evaluation. The results, however, do not change seriously. 
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From Fig. 7 - 8, we make a simple assumption that the symmetry of the joint pdf Pi2(^i, x 2 ) 
is represented as a regression line fitted by least-square method as follows 

logio x 2 = 9 log 10 xi + log 10 a . (4) 

In this form, the detailed balance (2) has the special symmetry, 9 = a = 1. For each scatter 
plot, we measure 6, a in the same range where Pareto index [i is estimated and the result is 
shown in Fig. 9. 

From this symmetry (ax\ <-> x 2 ), we extend the detailed balance (2) to 

fx 

Pi2(x 1 ,x 2 ) = Pi2( — ,ax! e ). (5) 



a 



We call this law the detailed quasi-balance. 

3 Pareto's law with varying Pareto index 

In this section, we derive Pareto's law with varying Pareto index by using the detailed quasi- 
balance (5). In the proof, we assume Gibrat's law (3) in the large scale region, because the 
number of data points is insufficient to observe Gibrat's law. 

Due to the relation of Pi 2 (xi, x 2 )dx\dx 2 = Pir(x±, R)dx\dR under the change of variables 
from (xi,x 2 ) to (x±,R), these two joint pdfs are related to each other 

P 1R (x 1 ,R) = x 1 e P 12 (x 1 ,x 2 ) , (6) 

where we use a modified ratio R = x 2 jx\ . From this relation, the detailed quasi-balance (5) is 
rewritten in terms of R) as follows: 

P 1R (x 1 ,R) = aR- 1 P 1R ((-j) ,a 2 R~ l ). (7) 

Substituting the joint pdf Pir(xi,R) for the conditional probability Q(R\xi), the detailed quasi- 
balance is expressed as 

Pi(xi) aQ{a 2 R- 1 \{x 2 /a) l l e ) 

1°) 



Pi{(x 2 /a) l ' e ) R Q(R\xi) 

By expanding Eq. (9) around R = a, the following differential equation is obtained 

a G'(a) 9 Pi(xi) + xi P'(zi) = . (10) 

The solution is given by 

P 1 {x l )=C l x l " aG '^ e . (11) 
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Here we consider two cumulative probability distributions P\(> X\) oc x\ m and P 2 {> 22) oc 
X2~^ 2 -, which lead 

Pi (ari) = Ci xf 1 " 1 , (12) 

P 2 ( X2 ) = C 2 x^- 1 . (13) 

From Eq. (11), (12) and (13), the relation between fi 2 and 9 is expressed as 

^±1 = 9. (14) 

This is an equation between detailed quasi-balance and Pareto's law in the non-equilibrium 
dynamical system. We confirm that the empirical data satisfy this correlation in Fig. 10. 



4 Conclusion 

In this study, we have investigated the dynamical behavior of non-equilibrium system in the 
large scale region by employing data on the assessed value of land in 1983 - 2006 Japan. We 
have identified the detailed quasi-balance (5) in the database, and have derived Pareto's law with 
varying Pareto index by assuming Gibrat's law (3). As a result, we have obtained a relation 
between the change of Pareto index \i and the parameter 9 in the detailed quasi-balance. The 
relation (14) has been confirmed in the empirical data nicely. 

What does the other parameter a mean? Because we demand detailed quasi-balance in the 
system, the area above the regression line (4) equals the area below it. The two parameters 9 
and a are, therefore, related to each other. The relation is expressed as 

= l-^log lo a. (15) 

Here 10 r is sufficient large number compared with the upper bound (10 7 ) where 9 and a are 
estimated. This is the reason why the two parameters 9 and a vary in opposite direction in 
Fig. 9. The relation (15) is confirmed directly in Fig. 11 where we set V to be 10. Consequently, 
the detailed quasi-balance has one parameter in principle. 

We should comment on several separations between 9 and (/xi + l)/(^2 + 1) in Fig. 10. An 
abrupt jump of Pareto index between 1984 and 1986 (2001 and 2002) is observed in Fig. 6. This 
means that the system changes vigorously in this period, where the symmetry is not represented 
as the regression line (4). Nevertheless, the dynamical equation (14) is valid in almost all the 
other quasistatic periods. 

For the next step, we should investigate the dynamical behavior in the middle scale region. 
In Ref. [8], it is reported that Pareto index jjl changes annually whereas the distributions in the 
middle region are stationary in time by analyzing income data for 1983 - 2001 USA. This phe- 
nomenon is explained by the Fokker-Planck equation [9] under two assumptions with respect to 
the change of income. The middle scale distributions of land prices, however, are not stationary. 
The distributions do not collapse onto a single curve by the normalization of the average land 
price. This difference is thought to be caused by the difference between the trend of increasing 
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(decreasing) income and the trend of increasing (decreasing) land price. In order to study the 
dynamical behavior in the middle scale region, we must identify each peculiar feature of middle 
scale distributions in the database [10]. 
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Figure 1: Cumulative probability distribution P(> x) of land prices in 1983 - 1985. The number 
of the data points is "16,975" in all cases. Data points are equally spaced in logarithm of land 
price. 
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Figure 2: Cumulative probability distribution P{> x) of land prices in 1986 - 1991. The number 
of the data points is "16,635", "16,635", "16,820", "16,840", "16,865" and "16,892", respectively. 
Data points are equally spaced in logarithm of land price. 
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Figure 3: Cumulative probability distribution P(> x) of land prices in 1992 - 1996. The number 
of the data points is "17,115", "20,555", "26,000", "30,000" and "30,000", respectively. Data 
points are equally spaced in logarithm of land price. 
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Figure 4: Cumulative probability distribution P{> x) of land prices in 1997 - 2001. The number 
of the data points is "30,300", "30,600", "30,800", "31,000" and "31,000", respectively. Data 
points are equally spaced in logarithm of land price. 
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Figure 5: Cumulative probability distribution P(> x) of land prices in 2002 - 2006. The number 
of the data points is "31,520", "31,866", "31,866", "31,230" and "31,230", respectively. Data 
points are equally spaced in logarithm of land price. 
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Figure 6: Annual change of Pareto index jj, from 1983 to 2006. For each year, Pareto index fi 
is estimated in the range of land prices from 2 x 10 5 to 10 7 yen/m 2 by using least-squares fit 
to data points equally spaced in logarithm of land price in cumulative probability distributions 
(Fig. 1-5). 
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Figure 7: The scatter plot of all pieces of land assessed in 1995 (xi) and 1996 (2:2). The number 
of data points is "11,278". 
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Figure 8: The scatter plot of all pieces of land assessed in 2002 (xi) and 2003 (22). The number 
of data points is "6,839" . 
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Figure 9: Annual change of 9 and a of Eq. (9) in the year (x 1 ,x 2 ) = (1983, 1984) - (2005, 2006). 
Because we demand the detailed quasi-balance (5) in the system, the two parameters 9 and a 
change in opposite direction. 
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Figure 10: Annual change of 9 and (jx\ + l)/(/X2 + 1) in the year (xi,X2) = (1983,1984) 
(2005,2006). 
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Figure 11: Annual change of 6 and 1 — 2/riog 10 a in the year (x\,X2) = (1983,1984) - 
(2005,2006). 
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